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Tagged jet measurements in high energy hadronic and nuclear reactions provide constraints on 
the energy and parton flavor origin of the parton shower that recoils against the tagging particle. 
Such additional insight can be especially beneficial in illuminating the mechanisms of heavy flavor 
production in proton-proton collisions at the LHC and their modification in the heavy ion envi¬ 
ronment, which are not fully understood. With this motivation, we present theoretical results for 
isolated-photon-tagged and B-meson-tagged b-jet production at ^/sNN = 5.1 TeV for comparison 
to the upcoming lead-lead data. We find that photon-tagged b-jets exhibit smaller momentum im¬ 
balance shift in nuclear matter, and correspondingly smaller energy loss, than photon-tagged light 
flavor jets. Our results show that B-meson tagging is most effective in ensuring that the dominant 
fraction of recoiling jets originate from prompt b-quarks. Interestingly, in this channel the large 
suppression of the cross section is not accompanied by a significant momentum imbalance shift. 


I. INTRODUCTION 

High transverse momentum (pr) production of energetic particles [l| and jets 0 in proton-proton (p-bp), proton- 
nucleus (p+A) and nucleus-nucleus (A-l-A) collisions plays an essential role in understanding the fundamental theory 
of strong interactions, Quantum Chromodynamics (QCD), and in probing the properties of cold nuclear matter 
and the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC) [1,0|- Recent theoretical and phenomenological advances, ranging from 
improvements in energy loss models @ to a unified treatment of vacuum and medium-induced parton showers 
have resulted in a good description of the light hadron cross sections attenuation in heavy ion collisions relative to 
binary collision-scaled p+p result, a phenomenon known as jet quenching. Predictions for light flavor reconstructed 
jet suppression Q were also confirmed by LHC measurements Q. However, early expectations that the corresponding 
suppression of D-meson and B-meson production will be significantly smaller due to heavy quark mass and color charge 
effects 0 were not supported by RHIC data EEl- The observed discrepancy between theory and experiment has 
inspired the development of a variety of theoretical models [IM3- Consistent understanding of light and heavy flavor 
suppression still remains a challenge for theory in light of the LHC results [l^ . 

Heavy flavor jets at the LHC have emerged as a new tool to further test the theory of heavy quark production, 
parton shower formation initiated by prompt b-quarks, and the modification of these processes in nuclear matter due 
to energy loss effects. Detailed studies of b-jets may also become part of the RHIC experimental program in the 
future |l^. In a previous work we studied single inclusive b-jet production in heavy ion collisions. The predicted 
suppression of the b-jet cross section in lead-lead (Pb-bPb) reactions is consistent with subsequent experimental 
measurements by the CMS Collaboration and supports a coherent energy loss picture, where heavy quark mass 
effects disappear at high px- At the same time, the inclusive nature of this observable does not allow to extract in a 
model-independent way the fraction of the b-jet energy that is dissipated in the medium. Furthermore, inclusive b-jets 
receive a large contribution from prompt light quarks and gluons, where heavy flavor emerges from parton splitting 
only in the late stages of the parton shower evolution. Thus, the connection between b-jets suppression and b-quark 
energy loss can be quite indirect. 

In this paper, our main goal is to identify differential observables that enhance the fraction of b-jets that originate 
from prompt b-quarks and provide constraints on the energy of the parton shower. To this end, we extend our 
calculations of inclusive b-jets to isolated-photon-tagged and B-meson-tagged b-jets production in heavy ion collisions 
at center-of-mass energy per nucleon-nucleon pair ^ysNN = 5.1 TeV, to be compared with experimental data from 
the forthcoming Pb-|-Pb run at the LHC. We discuss the advantages of these final states for studying the physics of 
b-quark production and propagation in dense QCD matter and present quantitative predictions for the tagged b-jet 
cross section suppression and the modification of the related momentum imbalance distribution in the QGP. We note, 
that earlier studies of photon-tagged [2^ - 1^ and heavy meson-tagged [26l - [^ heavy flavor in p-|-p and A-l-A reactions 
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FIG. 1. Left: The 7 +b-jet differential cross sections in p+p collisions at the Tevatron ^/s = 1.96 TeV is shown as a function 
of in two photon rapidity regions, \y'^\ < 1.0 (filled circles) and 1.5 < |y^| < 2.5 (open circles, multiplied by 0.3 for 
presentation), as measured by the DO Collaboration [^. The b-jet kinematics cuts are \y^\ < 1.5 and pip > 15 GeV. Red 
curves are from Pythia 8 simulations, using the anti-fcr algorithm with R = 0.5. Right: The same differential cross section 
from Pythia 8 simulation in p+p collisions at LHC yfs = 5.1 TeV as a function of Pp. The tagging particles (either the leading 
B-meson or the isolated 7 ) have the following kinematic cuts < 1.5 and Pp > 30 GeV. The b-jets are constructed using 

anti-lcT algorithm with R = 0.3 and have rapidity \y^\ <2. The blue (red) curve is for the B-meson-tagged (isolated- 7 -tagged) 
b-jet production. 


have also advocated the use of more exclusive final states to unravel the dynamics of heavy flavor production in 
hadronic and nuclear reactions. 

The rest of our paper is organized as follows. In Sec. II we present the evaluation of the isolated-photon-tagged 
and B-mesoir-tagged b-jet cross sections in p+p collisions using Pythia 8 simulations. In Sec. Ill we first describe 
the simulation of the medium-induced parton shower that recoils against the taggiirg particle and the related heavy 
quark energy loss effects. We then preseirt our phenomenological results for the tagged b-jet productioir in Pb+Pb 
collisions at the LHC. We conclude our paper in Sec. IV. 


II. PHOTON-TAGGED AND B-MESON-TAGGED B-JET PRODUCTION IN P+P COLLISIONS 

In this section, we present the evaluation of the differential cross sections for photon-tagged and B-meson-tagged 
b-jet production in p+p collisions using Pythia 8 [ 2 ^. This generator utilizes leading-order perturbative QCD matrix 
elements for 2 —>■ 2 processes, combined with a leading-logarithmic pT-ordered parton shower, and the Lund string 
model for hadronization. Specifically, we use the CTEQ6L1 parton distribution functions [s^ in the simulation and 
perform the jet clustering with the anti-Zcy algorithm j31| and a jet radius parameter R. In the photon-tagged b-jet 
production simulation, we hrst select the photon according to the desired kinematics and isolation cut specified below, 
and then we identify the b-jet on the away side with \(j)p — 4)j\ > 3 / 47 r, where 4>p {(j)j) is the azimuthal angle of the 
photon (b-jet). The b-jets are defined to be jets that contain at least one b-quark (or b-quark) inside the jet cone: 
a b-quark (or b-quark) is assigned to a jet if the radial separation from the reconstructed jet axis satisfies AR < R, 
where AR = y/(A^)^ + (Ay)^ with </> and y being the azimuthal angle and rapidity, respectively. For the B-meson- 
tagged b-jet production, which has a much larger cross section in comparison to 7 -tagged b-jet production, we select 
only the leading B-meson in each event. 

Measurements of the production cross section of isolated photons associated with b-jets in proton-antiproton (p+p) 
collisions at ^/s = 1.96 TeV were performed by both DO [s^ and CDF [s^ Collaborations at the Tevatron. Our baseline 
simulation can, thus, be validated against the experimental data, as shown in Fig. [T] (left). Here, the differential cross 
section da/dpp is plotted as a function of photon transverse momentum pp. We implement the experimental cuts 
in our simulation: for the b-jets \y^\ < 1.5, pip > 15 GeV, and the jet radius parameter R = 0.5. For photons, 
two rapidity cuts are considered: \y'^\ < 1.0 and 1.5 < \y^\ < 2.5. The photon isolation cut is also imposed in the 
simulation, which requires the hadronic transverse energy in a cone of radius Riso = \/{y — 2 / 7 )^ + (^ ~ ^ 7 ) = 0.4 
around the photon direction to be less than 7% of the photon energy. As can be seen in Fig. [T] (left), our simulations 
give a good description of the DO experimental data. 

In the following section we will present results on the nuclear modification of both isolated-photon-tagged and 
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FIG. 2. The fractional contributions of different subprocesses to the photon+b-jet production cross section are plotted as a 
function of the b-jet transverse momentum (left) and the photon transverse momentum (right) in p+p collisions at 
= 5.1 TeV at the LHC. The b-jet is clustered using anti-ftr jet algorithm with R = 0.3. We have integrated over the b-jet 
(photon) rapidity \y^\ < 2 (|p^| < 1.5), and > 30 GeV. 


B-meson-tagged b-jet production in Pb-|-Pb collisions at = 5.1 TeV. We first discuss the baseline calculation, 
i.e. the differential cross sections in p-l-p collisions at the same center-of-mass energy. In this simulation (and all 
the calculations below) we integrate over the b-jet rapidity \y^\ < 2 and the isolated-photon and B-meson rapidity 
|y 7 ,S| ^ With the transverse momentum > 30 GeV, the differential cross sections are plotted as a function 
of the b-jet transverse momentum pip in Fig. [T] (right). As can be seen from the figure, the y-tagged b-jet production 
cross section is comparable to the one measured at the Tevatron, while B-meson-tagged b-jet cross section can be more 
than two orders of magnitude larger. Both of those observables should be experimentally accessible at the LHC [sj. 
We note that an uncertainty in the absolute cross sections, for example due to different Pythia tunes, will cancel in 
the nuclear modihcation which we study in the next section. To reduce the statistical uncertainty/fluctuations in the 
baseline p+p cross sections, we have simulated around 10® events for both B-meson-tagged b-jets and y-tagged b-jets. 

In order to study energy loss effects and the medium-induced parton shower in heavy ion collisions, we further 
need detailed information on the flavor origin of b-jets produced in the p+p collisions at y/s = 5.1 TeV. Heavy quark 
production in hadronic collisions arises from various elementary 2 —>■ 2 hard partonic scattering channels, which can 




FIG. 3. The fractional contributions of different subprocesses to the b-quark+b-jet production cross section are plotted as a 
function of b-jet transverse momentum pip (left) and b-quark transverse momentum pp (right) in p+p collisions at y/s = 5.1 
TeV at the LHG. The b-jet is clustered using anti-fer jet algorithm with R = 0.3. We have integrated over the b-jet (b-quark) 
rapidity \y^\ < 2 (|y^| < 1.5), and > 30 GeV. 
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be easily separated in the Pythia simulations. In isolated-photon+b-jet production, for example, we have bg —>■ 76 , 
qg 75 , and qq —>■ jg. In the first subprocess, the produced b-quark initiates the b-jet. In heavy ion collisions, the 
medium modification of such b-jet production has a direct connection to the physical heavy quark energy loss (mass 
mb). On the other hand, the latter two cases can also generate b-jets through the parton shower. For the b-jets 
initiated through qg —?► 75 or qq —>■ 75 , the medium modification would resemble that of a massive quark (color-triplet 
state) or that of a massive gluon (color-octet state) of effective mass 2mb, since a bb pair is produced in the shower. We 
will use subscript symbol “(s)” to represent these different final states in the next section. Let us denote their relative 
contributions as Rb, Rq, and Rg, respectively, and plot them as a function of pip (left) and pp (right) in Fig. [21 As 
one can see from the figure, the j + b contributes about 50% of the cross section when the jet (and photon) transverse 
momentum are not very large pif ^ 100 GeV, and the contribution Rb gradually decreases as these momenta increase. 
On the other hand, the contribution from ^ + g gradually increases while the 7 -I- 5 contribution remains constant as 
pif^ increases. 

In this study, we focus on energetic boosted B-mesons, where the fragmentation process occurs predominantly 
outside of the medium and the interactions in QGP that we will consider happen at the partonic level [13; 113 ■ For 
B-meson-tagged b-jet production we first simulate the b-quark-fb-jet cross section in Pythia 8 . Once we implement 
all the interactions with the dense medium, we apply the b ^ B meson fragmentation function to convert the b-quark 
to a B-meson. At the partonic level, b-quark-fb-jet production arises from: (1) qq — >■ bb, gg —>■ bh; (2) gg gg, 
qq gg; (3) qg — qg; (4) qq — >■ qq, gg qq, qq — >■ qq. For the first case (1), we have bh in the final states, thus one 
b-quark initiates the b-jet, while the other b-quark fragments into the B-meson. This subprocess represents, on both 
the near side and the away side, the energy loss of a b-quark. For the remaining cases, the b-quark is produced in the 
parton shower and, thus, will mimic the energy loss of a massive gluon or a massive quark (see above). We implement 
the energy loss for such states in the formalism accordingly. The relative contributions of the various channels are 
presented in Fig. |3|as a function of pip (left) and Pp (right), respectively. As one can see from this figure, the 6 - 1-6 
contributions dominate in the whole kinematic region, which implies that the B-meson-tagged b-jet production has a 
much more direct connection to prompt b-quark energy loss. 


III. MODIFICATION OF TAGGED B-JET PRODUGTION IN DENSE QCD MATTER 


In this section, we present the medium-modification of 7 -tagged and B-tagged b-jet production in the ambiance 
of the QGP. Such final-state hot medium effects include radiative energy loss, caused by medium-induced parton 
splitting, as well as the dissipation of the energy of the parton shower through collisional interactions in the strongly 
interacting plasma MM,. The medium-modified 7 +b-jet cross section per binary nucleon-nucleon scattering is 


calculated as follows 20 , Mil: 


1 

(Nhin) dppdpip 



( 1 ) 


where (A^bin) is the average number of binary nucleon-nucleon interactions, the subscript “(s)” represents different 
final partonic states as described in Sec. m and P(s)(e) is the probability distribution that a fraction e of the hard 
scattering parton energy is converted to a medium-induced parton shower for a system with color state “(s)”. The 
hard production points are distributed according to the binary nucleon-nucleon collision density and the jets propa¬ 
gate through the medium that follows the participant number density and undergoes Bjorken expansion, which are 
determined by an optical Glauber model with an inelastic nucleon-nucleon cross section Um = 70 mb at ^/s = 5.1 
TeV d^. On the other hand, J(s)(e) = is a phase space Jacobian 

J(,)(e;i?,a;“'i) = {l - [l - /(,)(i?,a;“'i)] e}”' , (2) 

where {R, 0 ;'^°^*) is the fraction of the medium-induced parton shower energy that is retained inside the jet cone 
of radius R, as opposed to “lost” outside. Such a fraction is evaluated from the medium-induced gluon distribution 
(w is the energy and r is the angle) [lO, |4l| : 
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FIG. 4. Left: The nuclear modification factor Raa for 7 +b-jet production in Pb+Pb collisions at LHC ^/s = 5.1 TeV is plotted 
as a function of the b-jet transverse momentum p!j, and the photon Rapidities are integrated over |y^| < 1.5 and \y^\ < 2. 
Right: Raa for B-meson-tagged b-jet production as functions of b-jet py and B-meson p®. The rapidity cuts for B-meson and 
b-jet are \y^\ < 1.5 and \y^\ < 2. The coupling between the jet and the medium is Pmed = 2.0, and the mass for the collimated 
propagating of parent parton system for the b-jet is set to be M = 2 mj,. In these figures the suppression {Raa < 1) and 
enhancement {Raa < 1) are color coded, which can be seen in the online color version. 


where 0 ;'^°^* simulates the effects of collisional energy loss discussed in 


. It is obtained as a solution to 


'0 


dr / duj ^ = y = AEl 

duidr ^ 

iGshower 


icoll 

shower' 


( 4 ) 


Here, is part of the energy of the medium-induced parton shower dissipated in the QGP through collisional 

processes that was obtained in . To generate the medium-modified B-meson-|-b-jet cross section per binary collision, 
we first calculate the medium-modified b-quark-|-b-jet cross section as follows [3, 
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da 


AA 


(fVbin) dpypij- 
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t — e ./n 
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( pt /(1 - P)^J{sy)PT) 


dp\,dplp 


( 5 ) 


where (s) and {s') are the final-state parton color states in the hard partonic scattering level in the Pythia simulation, 
that further transmute into the b-jet and b-quark (this b-quark later on hadronizes into the B-meson), respectively. 
We then use the b ^ B fragmentation function Db^B{z) to fragment into the desired B-meson-tagged b-jet case: 


da 


AA 


dz 


(A^bin) dp^dpi. 


Db^B{z) 


da 


AA 


(A^bin) dp%,dpip 


( 6 ) 


Here, the momentum fraction z = where Pp is the B-meson transverse momentum, and we use the B-meson 

fragmentation functions as in [171 l2fi| | . 

Let us now turn to the phenomenological results for tagged b-jet production in Pb-|-Pb collisions at the LHC. Our 
simulations are performed at a center-of-mass energy per nucleon-nucleon pair ^/smn = 5.1 TeV for direct comparison 
to the planned experimental measurements at the LHC. We focus on the most central collisions with average number 
of participants Apart = 360 and the central rapidity region with photon or B-meson rapidity < 1.5 and b- 

jet rapidity \y^\ < 2. We choose the b-jet reconstruction radius R = 0.3 and require back-to-back configuration 
107,3 ~ 0il > 3 / 47 r. We implement the isolation cut for the photon that requires the hadronic transverse energy in 
the cone of radius R\so = 0-4 around the photon direction to be less than 7% of the photon energy. 

The many-body QCD dynamics that modifies the tagged b-jet production in heavy ion collisions will be manifest 
in the deviation of the A-l-A cross section from the binary collision scaled p-bp results. 
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( 7 ) 


Theoretical predictions for the nuclear modification factor Raa a-re presented in Fig. 21 for 7 -l-b-jet production (left) 
and B-meson-l-b-jet production (right), respectively. We include both collisional and radiative energy loss effects, and 
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FIG. 5. The isolated-photon-tagged b-jet asymmetry distribution is shown for different coupling strengths between the jet and 
the medium. We have b-jet transverse momentum > 40 GeV. Left panel: 40 < < 80 GeV. Right Panel: 80 < p^ < 120 
GeV. Bands correspond to a range of masses of the propagating system between mb and 2mf,. Blue, red, and green bands 
correspond to (?med = 1 . 8 , 2 . 0 , 2 . 2 , respectively. 




FIG. 6 . The B-meson-tagged b-jet asymmetry distribution is shown for different coupling strengths between the jet and the 
medium. Kinematic cuts and jet-medium couplings are the same as in Fig. [S] with the replacement 7 —^-B. 


the mass M of the propagating parton system in the final state is implemented as in [ 2 ^. In this illustration for 
the 3D plots of i?^^we set M = 2mb, twice the b-quark mass, and choose the coupling between the jet and the 
medium t/med = 2.0 [1, The coupling between the final-state partons and the medium is the largest uncertainty 
in our nuclear modification model. It affects the absolute value of the cross section suppression/enhancement but 
does not alter the qualitative behavior of Raa, which we now discuss. As one can see clearly, the largest suppression 
is observed when the transverse momenta of the trigger particle and the b-jet are similar; ^ « p^. For p^ < p^, 
there can be an enhancement of the cross section for 7 -l-b-jet production, which is consistent with what was observed 
in 7 -|-light-flavor-jet production 1^. For p^ > p^, the b-jet suppression is somewhat smaller than the quenching 
of inclusive b-jets [^, [^, because the fraction of 7 -tagged b-jets originating from a prompt gluon is smaller. The 
fluctuations arise from the Pythia 8 limited statistics for small, steeply falling baseline cross sections, especially when 
and py are significantly different from each other. On the other hand, we find that the overall suppression is 
stronger in B-meson-|-b-jet production than that of photon-|-b-jet production since both the parton that fragments 
into the B-meson and the b-jet lose energy when they traverse the medium, while the isolated-photon does not lose 
energy. When p^ ^ P-f , or pip <C Py there can be a strong enhancement in Raa 0 , which we don’t see in the more 
limited kinematic range that we study. The overall 2D suppression pattern is more uniform than for light jets. 

Next, we consider the tagged b-jet event asymmetry. For this purpose, we define the following variables to represent 
the momentum imbalance of the 7 -l-b-jet as well as B-|-b-jet events, Zj^ = pip/pi^, ZjB = Pt/Pt- Based on our 
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TABLE I. Theoretical results for {zj-y) {{zjB)) in photon-tagged (B-tagged) b-jet production in p+p and central Pb-|-Pb 
reactions at center-of-mass energy ^sjvat = 5.1 TeV at LHC. For the photon (B-meson), we have two kinematic cuts, while 
pij, > 40 GeV for the b-jet. 


System 

(zjj) 

(zjb) 

40 < Pp < 80 GeV 

80<pI< 120 GeV 

40 < pI < 80 GeV 

80 < pf < 120 GeV 

P+P 

A+A, g = 1.8, M = 2mh 

0.98 

0.95 

0.78 

0.74 

1.04 

1.06 

0.82 

0.83 

A+A, g = 1.8, M = rtib 

0.95 

0.74 

1.07 

0.83 

A+A, g = 2.0, M = 2mh 

0.94 

0.73 

1.07 

0.83 

A+A, g = 2.0, M = mt, 

0.94 

0.73 

1.09 

0.83 

A+A, g = 2.2, M = 2nih 

0.93 

0.71 

1.08 

0.83 

A+A, g = 2.2, M = mi, 

0.93 

0.71 

1.10 

0.83 


simulations of the double differential cross sections da/ dpip 


above, we evaluate the ^-distributions as follows: 


da 

dzj'r 


dpi 


p!p da(p7p = plplzj^) 


rp ^ 

Z- 
31 


dppdplp 


da 

dZjB 


dpL 


p!p da(p^ =pl^/zjB) 
ZjB dppdptp 


( 8 ) 


where p^™™ and are specihed by the experimental cuts on the b-jet transverse momentum. In our predictions 

for the future Pb-|-Pb run at the LHC we choose the b-jet transverse momentum pp > 40 GeV. On the other hand, 
we explore two different transverse momentum cuts for the 7 and the B-meson: (1) 40 < Pp^ < 80 GeV, (2) 
80 < Pp^ < 120 GeV. We see that cut (2) leads to more asymmetric distribution of 7 (B)-|-b-jet production than cut 
( 1 ). 

The momentum imbalance distributions are given in Fig. [5] for 7 -l-b-jet production, and in Fig. [ 6 ] for B+b-jet 
production, respectively. Note that here we present the absolute differential cross section da/dz instead of the usually 
measured normalized asymmetry distribution {lla)daldz 0,111, SMI. In both plots the left panel corresponds 
to cut selection (1), while the right panel corresponds to cut (2) - the more asymmetric one. The black histogram 
corresponds to the p+p baseline. The colored histograms correspond to the cross sections per binary collision in 
central 0-10% Pb+Pb collisions, with bands representing a range of masses of the propagating system between mb 
and 2mb. Blue, red, and green bands correspond to pmed = 1-8, 2.0, 2.2, respectively. The range of these values has 
worked well in describing and/or medicting the dijet 0 and photon-tagged jet [111 asymmetry distribution, as well 
as inclusive b-jet [ 2 ^ and hadron [1] suppression at the LHG. 

For the 7 +b-jet case in Fig. [5] and the more symmetric cut (1) we see a moderate suppression of the cross section 
which is strongest when pp ~ pp and zj-y ^1. At the same time, the final-state interactions in the dense medium 

will shift the asymmetry variable Zjy down to smaller values. This arises from the energy loss of the b-jet, while 

the photon escapes out of the medium unscathed. On the other hand, for the more asymmetric cut (2), we find an 
enhancement for Zjy ^ 0.6 and a change in the shape of the imbalance distribution. These features are all consistent 
with the asymmetric suppression pattern seen in the 3D plot of Raa in Fig.|Tl Note that the fluctuation/uncertainty 
in Raa should only have minimal influence in our plot here since they happen mainly in the region when pp and pip 
are significantly different from each other (by a factor of 3 or more), which corresponds to blow 1/3 or above 3 
and has much smaller cross sections. To further quantify such down-shift in Zjy distribution [jj, we define the mean 
value of Zjp (similarly also for Zjb) as 

{^31) = (/ / (/ = (/ / (/ 

and show its value (zjy) in Table. U for 7 +b-jet production. From these calculations, we see there is around 5% (10%) 

down-shift in (zjy) for the more symmetric cut ( 1 ) (more asymmetric cut ( 2 )), which is somewhat smaller than those 
values (^ 15%) in the 7 +light-jet production, as studied in [s^ and measured in [4^ . 


Here and in the rest of the paper, down-shift or up-shift are always relative to the p+p baseline. 
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On the other hand, for the B+b-jet case in Fig. |6]we find strong suppression for both cuts, (1) and (2). For the 
more symmetric cut (1), as can be seen in Fig. [5] (left), the distribution is peaked slightly above ZjB = 1- This is 
also illustrated by the average value {zjb) shown in Table HI The reason for such behavior is that it is the py of the 
b-quark (bigger than the of the B-meson) that balances the transverse momentum of the away-side jet. Besides 
suppression, in contrast to y+b-jet production there is no noticeable down-shift of the asymmetry ZjB distribution. 
The calculated (zjs) for both cuts even show a slight increase in going from p+p to Pb-|-Pb collisions, as seen from 
Table HI This can be understood as follows: as we discussed in the previous section, the b + h final state dominates 
b-|-b-jet production, recall Fig. |31 If we consider this dominant b-\-h channel and inspect Eq. ([5]), on average the b-jet 
loses less energy, (1 — /(s)(i?, w'^°"))e, than the b quark, e. This, in turn, leads to a slight up-shift of ZjB to larger 
values. This behavior is very different from the light dijet asymmetry enhancement observed in experiments (^. [d^. 
which corresponds to a momentum imbalance down-shift. It suggests that the large asymmetry enhancement for light 
dijets is likely driven by the different energy loss for quarks and gluons, proportional io Cp and Ca, in the dominant 
q-|-g hnal-state channel. The value of the up-shift reported in Table HI might have a small numerical dependence on the 
details of fragmentation into B-mesons, but this will not change the qualitative difference in the behavior of {zj-y) and 
{zjb)- It will be important in the future to use tagged b-jet Pb+Pb measurements at the LHC to elucidate the precise 
physics mechanisms responsible for the large nuclear modification of heavy flavor production in heavy ion collisions. 


IV. CONCLUSIONS 

In summary, we presented theoretical predictions for the 2D nuclear modihcation Raa and the related momentum 
imbalance shift of isolated-photon-tagged and B-meson-tagged b-jets in Pb-|-Pb collisions at ^/snn = 5.1 TeV at 
the LHC. We validated the Pythia 8 simulations of tagged b-jet cross sections in nucleon-nucleon collisions through 
comparison of photon-tagged b-jet results with p+p data at -/s = 1.96 TeV from the Tevatron. We also found that 
particle tagging can significantly increase the fraction of recoiling b-jets that originate from prompt b-quark relative to 
the inclusive b-jet case. While in the latter case this fraction can be as low as 20%, B-meson tagging in particular can 
increase the contribution of prompt b-quarks to 70-80%. In Pb+Pb collisions in the LHC we further considered the 
medium-induced parton shower in the soft gluon energy loss limit and any additional dissipation of its energy in the 
QGP due to collisional interactions. We found significant nuclear suppression when the trigger particle momentum 
is similar to the b-jet momentum p^^ ^ pip. A modest 5 — 10% increase in the transverse momentum imbalance for 
7 +b-jet production in Pb+Pb collisions is predicted from our calculations, depending on the specific kinematic cuts, 
which is slightly smaller than that observed for y+light jet production. For B-tagged b-jets we found an even larger 
suppression of the double differential cross section since both the jet and the tagging b-quark lose energy. On the 
other hand, the asymmetry variable ZjB in B+b-jet production shows a slight increase from p+p to A+A collisions, 
a behavior quite different from the one observed in dijet asymmetry distributions. By using the flexibility of b-jet 
tagging, comparison of upcoming experimental measurements to theoretical calculations, such as the ones presented 
here, can provide us with unique new insights into heavy flavor dynamics in the strongly-interacting plasma. In the 
future, we expect that the theoretical uncertainty of b-jet production in nuclear matter can be further reduced by 
going beyond the current radiative energy loss approximations im and by including higher order computations for 
multiple scattering [i^, with full account of heavy quark mass effects. 
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